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REMOTE SENSING OF LAND USE AND WATER QUALITY 
RELATIONSHIPS WISCONSIN SHORE. LAKE MICHIGAN 

by 

R,K. Haugen. H.L. McKim and T.L. Mariar 


INTRODUCTION 

The Great Lakes are one of the nation's inosi valuable resources and have been a major factor In 
the settlement and economic grovirth of the midcontinental area of both the United States and 
Canada. The lakes are utilized for many purposes, such as recreation, transportation, waste disposal, 
water supply, hydroelectric power producUon, and commercial fishing. A gradual deterioration in 
water quality caused by these uses has been documented (Beeton 1970), with Lakes Erie and 
Ontario showing the greatest amounts of pollution and eutrophication. Lake Michigan, the focus 
of the initial phase of this study, does not exhibit the advanced deterioration of the eastern lakes, 
but there are serious problem areas in Green Bay and at the southern end of the lake, as well as indi- 
cations of long-term deterioration in the open waters (Beeton 1970). The seriousness of this con- 
dition in Lake Michigan is compounded by iu low water turnover rate, compared to Ote eastern 
lakes, and most of the major streams which feed Lake Michigan are already polluted. 

As lar d use patterns evolve, generally following a progression from forest to agriculture to urban 
or industiial uses, stream flow and water quality characteristics also change. The U.S. Army Corps 
of Engineers, along with its responsibility for the maintenance of Inland navigable waterways, is 
concerned with the environmental impact of pollutants, especially as they affect lakes and streams. 
Stream sedimentation, a form of pollution which can often be related to economic activity, is of 
particular interest in this invesUgation. 

Objectives 

This investigation examines the feasibility of using remote sensing methods to assess, rai^dly and 
economically on a region^ scale, the effect of land use as it influences sedintent loading in streams. 

To accomplidi this objective, the following tasks were completed: 1 ) the acquisition and evalu- 
ation of a large variety of remote sensing imagery, 2) the development of techniques for mapping 
andquantifwation of land use patterns, 3) the anessmehtof the cost effectiveness of these mapping 
and measurement techniques, and 4) the collection and analysis of water quality samples for com- 
parison to land use charactMistics and correlation with signatures identifiable by remote smising 
tedmiques. 

Approach and projea history 

This study was initiated in April 1972. A field trip was uken in May 1972 for an overview of 
the test are.i and for the acquisition of data from sute and federal agencies in Wisconsin. The initial 


phdM! ol the »tudy was primarily devoted to acquiring published data and evaluating a variety of 
imagery provided by NASA in September 1972, 

NASA provided imagery for the project by utilizing two aircraft: the General Dynamics RB>57F 
(later designated the WB-S7I ), a mid*wing, four-jet engine airaaft capable of altitudes in excess 
of 60,000 It, and the Lockheed NP-3, a low-wing, four-engine aircraft designed to carry a variety of 
sensors at altitudes under 20,000 ft. Some imagery was also acquired by the investigators in a rented 
Cessna 1 85, a high-wing, single engine aircraft modified to accommodate ur aerial camera, Consid- 
er.ible effort was devoted to experimentation with mapping techniques by using two types of 
densitometry, fhe intent was to develop methods for the measurement of land use types directly 
from imagery. Using color infrared emulsions and filtering, some success was achieved in differenti- 
ating U.S, Geological Survey (USGS) Land Use Classification level I, but only where soil color was 
light. A system of densiiometric measurement with hand-drawn overlay maps was eventually de- 
veloped. This technique provided the necessary data for subsequent comparisons with water quality 
parameters. 

Two methods were used to extract from imagery water quality information that could be cor- 
related to held measurements. The first method was the development of a color classification for 
streams in the study area, based on a dcnsttomctric method. The purpose of the classification was 
to quantify photographic stream color differences displayed in the 1972 high altitude Ektachrome 
(SO-397) imagery. The results of this approach, however, were not sufficiently consistent to be of 
practical value. During the 1974 missions, ground truth data on reflected and incident radiation in 
the CRTS (Earth Resources Technology Satellite) MSS spectral bands were obtained. This method 
provided more reliable results and permitted significant correlation with measured concurrent sus- 
pended sediment toads for several major stf’cams in the test area. These data were then compared 
to measurements of land use and soil types in the individual watersheds for analysis of regional 
jrefationships. 

In July 1973, a chapter was prepared on land use mapping for the Corps of Engineers Handbook, 
Remote Sensing for Cnvlronment/il Analysis (Haugen and Splett 1974). This chapter, based on land 
use mapping procedures developed during dtc project, provides practical discussion and guidelines 
for land use mapping. ERTS-1 and high altitude color infrared imagery are compared for land use 
mapping applications within misuons of the Corps of Engineers, 


STUDY AREA 

The Wisconsin shoreline of Lake Michigan, NASA Test Site 31 1 (Fig. 1 ), was selected for this in- 
vestigation for several reasons: 1 ) the area is predominantly agricultural, widi relatively isolated con- 
centrations of urban and industrial land use, so that the contributions to pollution levels of the basic 
lutd use types can be readily assessed; 2) the general pollution of Lake Michigan is still low enough 
so that future benefits from identification of present and potential ‘‘trouble spots” can be realized; 

3) the streams entering Lake Michigan from the Wisconsin shoreline arc relati^'ely short, so that 
entire watershed areas can be mapped within about 50 miles of the shoreline. 

Physical setting 

The test area lies within the Eastern Ridges and Lowlands physiographic subprovince of Wiscon- 
sin. The eastern portion, the Eastern Ridges, is entirely underlain by the Niagara Cuesta limestone 
which forms the Door County Peninsula and extends to the southern boundary of Wisconsin, slop- 
ing gently toward Lake Michigan. The other major physiographic subprovince is the Lake Winnebago- 
Green Bay Lowland, a gently sloping plain which contains Lake Winnebago, the Fox River Valley, 



and a submerged portion forming Green Bay. The entire landscape of the test ata has been exten* 
sivety modified by Pleistocene glaciation* with glaciai deposits forming most of the local relief 
which ranges up to 200 ft. The Lake Michigan shoreline is formed from tills and lacustrine sedi* 
ments of higher glacial iakes^ resulting in wave^ut cliffs 100-^I2U ff high in much of the test area. 
Erosion of these unconsolidated materials continues at the present time, resulting in a large sedi* 
ment load in the nearshore waters. 

The major rivers of the study area are the Fox, Wolf, Menominee, and Oconto which drain into 
Green Bay, and the East and West Twin, Manitowoc and Sheboygan which drain into Lake Michi^n 
(from north to south along the shoreline). Hie rivers flowing directly into Lake Michigan do not 
have sources west of the Niagara Cucsta, with the single exception of the Manitowoc, which is 


believed tu have drained Lake Winnebago during glacial iimes« Many lakes dot the glacially modi* 
fied landscape, few of them with an a^ea exceeding 3000 acres. Lake Winnebago, 215 square miles 
in ar^a, is the largest lake in the state and was formed fce^’ind a morainal dam during glacial retreat. 

Sod ossfitraiions in the study area include the Niirthcrn and Lastern Sandy and Loamy Reddish 
Drift Uplands and i'lains, Northern Silly Uplands and Plains, Northern Loamy Uplands and Plains, 
Northern Sandy Uplands and Plains, and the Stream Bottoms and Major Wetlands in the area between 
the Oconto and Fox Rivers. The major soil associations in tltc Last and West Twin and Manitowoc 
River Watersheds are the Soils of tile Northern and Lastem Clayey and Loamy Reddish Drift Up- 
lands and Plains. Minor soil associations include well to poorly drained sandy soils and poorly 
drained dcpres.*^onal soils with some peats and mucks. 

The climate of the test area is continental, with long and severe winters, especially in the north, 
and very warm summers. Lake Michigan modifies temperatures in a narrow bell bordering the lake, 
but the effect is minimized by the genera) easterly and r juth-casterly trend of most weather system 
movements through the areas. The warming effects o"' Lake Michigan arc, therefore, much less for 
eastern Wisconsin than for western Michigan. 

January is the coldest month, with an average temperature ranging from 13®F in the north to 
24""F in the south. During July, the warmest month, average temperatures range from 66° to 73*F 
north to south. The fro5t4rce season ranges from 10 June to 30 August in the north, to 25 April 
to 20 October in the south. The wettest months arc May through September, and annual precipita- 
tion touls range from 29 to 33 in., north to south. Snowfall ranges from 60 in. in the north to 30 
in. in the south. 

Weather conditions in eastern Wisconsin prior to and during aircraft missions and water sampling 
operation ^ arc Indicated by the temperature and precipitation data in Appendix A. Generally, the 
high altitude RB-57F flights required fair weather, with few clouds; however, scattered rain showers 
preceded the flights in much of the lest area* The June 1974 NP-3 nighttime thermal scan was not 
completed due to heavy ground fog, and weather conditions during September and October 1973 re- 
sulted in two canceled missions, which were rescheduled and accomplished during May and June 
1974. 

Cultural setting 

The population within the test area is in the neighborhood of 3,000,000. The largest population 
concentration Is in and around Milwaukee, which is also the economic center of the area. Manufac- 
turing is by far the largest employer within the test area, followed by agriculture, forestry, consuuc- 
lion and mining. Over half the northern part of the test area is forest, and wood processing firms 
are a major part of the region's economy. Much of the wood processing industry is centered in the 
Fox River Valley. Mining activity primarily Includes extraction of clay, limestone, sand and gravel, 
with some iron ore and granite being mined In the north. The recreation industry is also important, 
particularly in the northern forested areas. 

Agriculture Is important throughout the area. Grain crops are dominant In the south, and vege- 
tables, small fruits and horticultural specialities arc raised in the central and northeast portions. 

Some irrigation is practiced with specialty crops. Dairying and related processing industries are 
found throughout the area. 
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DOCUMENTED POLLUTION SOURCES IN THE STUDY AREA 


There has been considerable atiiviiy in the identification and control of aquatic pollution in 
Citsiern Wisconsin. Major sources of pollution have been idcmilicdi and recommendations for re- 
medial action have been made by the Wisconsin Department of Natural Resources. Three divisions 
of the U.5* Department of Agriculture (the Soil Conservation Service, the Exonomic Research Ser- 
vice, and the U.S. Forest Service), arc engaged in a cooperative study of flooding and water manage- 
ment problems, the Southeast Wisconsin River Basin Survey. 

A map of known aquatic pollution sources (Fig. 2) was compiled from information gathered by 
the above agencies, especially the Wisconsin Department of Natural Resources, Environmental Pro- 
tection Division. The point sources of pollution arc mapped in five categories: 1) sewage, 2) dairy, 
3) cannery, 4) thermal and 5) other types. (As will be discussed in subsequent sections, few of these 
point sources arc visible in the small-scale photography acquired during Mission 205 or Mission 273.) 

The largest contributors to polluted waters arc inadequate sewage treatment systems. This prob- 
lem is not unique to the study area, but typifies other areas throughout the United States. Most «>f 
the major cities in the study area have sewage treatment plants, but many plants are working at 
greater capacities than they were designed to handle. Some smaller towns have no treatment plants, 
and inadequate septic tanks contribute to the problem. Attempts arc being made to reduce the ex- 
cessive amounts of pollutants reaching Uic streams. 

The dairy industry is one of the Icadin^^ sources of income within the study area. As a result, 
many dairy processing plants operate in the area, and in most eases, their wastes reach streams with 
little or no ireaimcnU Dairy farms also add to the pollution of the area because of wastes from 
their livestock and fields. In other instances, dairy wastes enter existing sewage treatment plants, 
contributing to the overloading of some of them. 

Vegetable and fruit production is an important occupation among farmers in the area; conse- 
quently there are canneries that contribute wastes to streams. These wastes add various chemicals 
from washing operations and increase the organic load. 

A limited number of thermal pollution sources exist within the study area. These include electric 
power generating facilities and sewage treatment plants. The mafor sources of waste heat arc two 
nuclear power plants, the 527,000-kW plant at Kewaunee and the 597,000-kW plant at Point Beach. 
These plants arc located on Lake Michigan and use water directly from the lake for cooling. 

The "other*' category of point sources of pollution includes paper-processing factories, which 
contribute both chemical and organic wastes, industrial plating facilities and tanneries. 


EVALUATION OF IMAGERY 


Photographic imagery 

The flight lines for tnc RB-57F photographic missions during 1 972 arc indicated in Figure 3. 
During 1973, a low altitude thermal scan (Mission 235), was flown with coverage of the coastal 
portions of the test area along with the Fox River {Fig. 4). An attempted second RB-57F flight in 
September 1973 was aborted due to unfavorable weather conditions over the test area. The 1972 
RB-57F coverage was repeated in 1974 during Missions 273 (3-7 May 1974) and 284 (12 June 
1974). Separate flights of these missions were required to cover the test area because of poor 
weather conditions* A second coverage with the NP-3 aircraft was also obtained on 31 May 1974 



”P 0 40hm 

FigurtZ Documented oquMk pollution tources in the study afta. 



f igurt 3. Cenerafiied flight fines for RB-5 ?r photographic missions. 


(Mission 272), which included the Wisconsin coast, the Fox, Manitowoc and Sheboygan Rivers* The 
flight lines for the 1974 photographic and thermal scanning missions are essentially the s^ime as shown 
in Figures 3 and 4. 

The RB-S7F aircraft was ec|uipped with two 9**n. format RC-8 cameras (each with a 6-in. lens), a 
9-in. format Zeiss aerial camera (with a !2-in. lens), a 70-mm Kasselblad camera (with 40-mm and 
80-mm lenses), and an RS-7 thermal scanner filtered for the 10.2- to 12J-|im spectral region. Imagery 
from the various cameras consists of 9-in. Ektachromc Type 50-397, 9-in. and 70-mm Aerochrome 
Infrared Type 2443, and a Mack and while positive of the RS-7 thermal scanner image. Sequential 
photography was obtained with the two RC-S cameras timed for 60K overlap, and the Zeiss camera 
limed for 22% overlap. As the RB-57F flew at altitudes of 59,000 to 59,800 ft above ground level, 
the approximate scale of the 9*in. photography was 1 : 1 20,000 for the 64n. lens and 1 :60,000 for the 
12-in. lens.* Table I shows the imagery acquisition schedule for the duration of the project. 


* The exact scale was 1:118,000 for the 6-in. lens and 1:59,000 for the 12-in. lens. 
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Figuft 4. Flight lines for low^titude thermal scannirtg missions. 


The qiulity of the photography v«ied from good to excellent in the 9*in. format, and from poor to 
excellent in the TO^nm format. Evaluation of each type of imagery obtained for this study should 
be helpful for planning future misfiom. 

Aerochrome infrared film (9Hn. fornat) was exposed to two cameras (RC-8 and Zei») with 6-in. 
and 1 2-in. lenses. The 1 :120,000 scale colorinfrared photography was used for direct comparison 
with the non-infrared Ektachrome Type SO-397, which was exposed at the tame scale. The t ;60,000 
photo^aphy was used to verify infonrution derived from the smaller wale photogrerV. Photography 
was also attempted in 70-mm format with Aerochrome infrared film and with a multispactral 4-camcra 
system. A description of the various photographic fomwts follows. 

Aerochrome infrared film, type 244X 94n. format, 64n. lens, xale 1:120,000, On the 1 :1 20,000 
Aerochrome infrared photofraphy, land use patterns are readily distinguished but individual uniu are 
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Tabic L Imagery acquisition schedule for NASA and for 
investigators of this study* 


Date 

Aircraft 

Mission 

Altitude 

M 

Sensor 

Film type! 
wavelength 

Imagery acquisition by NASA 





‘1 June 1972 

RB-57F 

205 

60,000 

RC-8, 6-in. lens 

2443 





RC-8, 6-in. lens 

SO-397 





Zeiss, 1 2-in. lens 

2443 





Hassciblad, 40-mm lens 

2443 





RS-7 - thermal 

10-12 jum 

24 June 1973 

NP-3 

235 

3,500-4,000 

RS-14 - thermal scanner 

8-1 4 /Jm 





PM IS 

K-band 





PRT5 


3-7 May 1974 

RB-57F 

273 

60,000 

RC-8, 6-in. lens 

2443 





RC-8, 6-in. lens 

SO-397 





Zeiss, 12-in lens 

2443 





Hassciblad, 40-mm and 

2443 





80-mm lenses 






Hasselblad, 4-camcra 

2402 





multispectral 40- and 






80-mm lenses 


31 May 1974 

NP-3 

272 

2,500-4,000 

PMIS 

K-band 





RS-14 

8-1 4 pm 

12 June 1974 

RM-57F 

284 

60,000 

Same as Mission 273 


Imagery acquisition by investigators 




3 Oct 1973 

Cessna 

185 

5,000 

Hassciblad, I00-mm lens. 

PX, 2402 





4-camcra multispectral 

2443 






EKMS, 2448 

7 May 1974 

Cessna 

185 

5,000 

Zeiss, 6-in. lens 

2443 

nv 


not easily discerned* Water bodies suc!j as lakes, ponds, rivers and streams arc easily seen* Patterns 
are observed in larger water bodies which may be indicators of waterborne pollutants* A wide dis- 
parity is noted in the color rendition of water bodies, and it is obvious that factors other than pollu- 
tion arc contributing to tonal appearance. Some of these factors arc depth of water, surface rough- 
ness, sun angle, and water coloration. This 1:120,000 scale photography is useful in preparing land 
use maps of a general nature and is excellent for long-term monitoring of urban expansion. However, 
the detail available «it this scale is not adequate to make positive identification of all observed features. 
The photography obtained in the 9-in. format with the longer focal length lens (12-in.) was used 
wherever greater detail was needed in the study. 

Aerochrome InfraredF/lm, type 2443, format, 124n, lens, scale 1:60,000, The detail ob- 
served in this photography is excellent; center markings on highways may be seen, and individual 
automobiles arc seen on highways and in parking lots. Vegetation Is readily separated as to type 
and si/e. A peculiarity of the imagery is :hc three-dimensional appearance of the red vegetation 
signature, a feature that is very helpful in determining relative vegetation heights. The effect is ap- 
parently caused by the diffraction characteristics of the red color in the viewing optics. Surface 
water is contrasted sharply with the background, and drainageways can be traced with relative ease. 
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1\uiems ippciUmi; in Krnuiul features lead to the spcLUlatiun (not investigated in tlte study) that 
relative surface moistiiu* differences in soils may he seen. Water depth is difficult to estimate hut 
is indicated hs changes in the blue rendition of the water* f r)r eJcamplc* a large swimming pool is 
sluiwn in the imagery with water dial is a very light hluc at one end and gradually liecnming dark 
blue al the opposite eiuL This color gradation isallected by other factors than depth, liowever, and 
is not a reliable inJkator in al) instances. 

U rndnomc fnfrarvJ, typv 7()-mm fnrmaf (7^ i-ni square}, uO- and 80-tvm femc% 12 + 30B 
tiUtr, this photography was taken on two days* On ^ Mav it was exposed with 40-mm lenses 
un the I lassellilad cameras and on 7 May 1974 it was exposed with SO*mm lenses on llic same cameras, 
f he usefulness of »he two days* phohigraphy varied t onsiderablv . I arge palicrns such as urban areas, 

\ age fakes nr ponds and streams could be secn» hut finer detail needed for land use mapping was 
missing. The probable cause ol such poor resolution was aimosplieric altcmialion coupled with 
exposure variances* The TOrntm formal was mure difficult to work with than larger formats and 
was considered the least desirable for land use mapping* 

liassclhfad andn-spec tra! imaqcry. The NASA d-camcra T lasselblad system was used w ilh 
Uralten 24 and 4S filters and 2402 and 2442 (IR) film Xn obtain black and white pholograpliy of 
the test areas which would approximate the various spectral hands obtained in the [' RTS*1 imagery 
(blue, green, red and lUMr IR)* The 2402 film with a Wratten 24 filter produced the best result lor 
single band viewing and was the only hand considered suitable for gross land use mapping (level 1 )* 

C olor reconstruction of selected scenes was attempted using the muUhsptclral viewer, but the use* 
fulness of this technique (optical mixing) was limited due to the quality ol the photography. It is 
possible to ernphasi/e natural and cultural features by using the multi-speclrai viewer, hut in no 
instance was more information obtained from this technique than could be obtained from the larger 
color transparencies. 

CoruJuvom^ It was found in this study that the aerial photography dues not give enough infor* 
nutiejn to allow direct monitoring of point sources of water pollution. It does, however, provide 
information on gross pollution in water bodies, primarily in those cases where the pollution con* 
trasls with the normal coloration of the water. The ability to do detailed land use and drain»ige 
mapping from this scale photograpiiy is its principal asset. Smaller scale photography is also useful 
in land use mapping, when supplemented with adequate larger scale photography for detail checking. 

Thermal imagery 

Patterns indicative of thermal pollution are generally recognized as thermal plumes emanating 
from point sources such as industrial areas, thermoelectric generating plants (with both nuclear and 
conventional fuels), and sewage treatment plants. The detection of thermal plumes is possible on 
small-scale thermal imagery (such as from an RS*7 scanner), but no quantitative assessment of actual 
lenipcraiurc differences is possible, since calibrated temperature data are not available. The evidence 
of thermal pollution is therefore based entirely on relative lone signatures; for example, the rendition 
of dark tones on a positive film indicates cool icmpcraiurus, while light tones indicate warm 
temperatures. 

RS-7 thermal imagery, Tlic iliermal imagery for this study was first obtained with the RS-7 infra- 
red scanner, wliich is sensitive to the 10*2- to 12*5-pm spectra! range, the far-infrared region of the 
spectrum. Most thermal imagery is acquired at night to avoid the effect of reflected radiation and at 
altitudes <inder 10,000 ft so that the resolution will be sufficient for the detection of thermal point 
sources. For this study, however, high-altitudc, daytime RS-7 imagery was obtained on an experi- 
mcnti^i basis. In terms of the obiectivcs of this study, liule evidence of thermal pollution patterns 
could be seen. The only clearly identifiable man-made plume on the original imagery was the plume 
of the Point Beach Nuclear Poiver Plant (located at point 4, F!g. 5) where thermal differences could 
be seen extending about two miles southward of the plant. 
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The* pUintes III scvml rivers ctueriii)* lake Miililiian amJ (ireen tSay incluilo the Milwaukee, the 
[ i»\, and the Menominee Kivers, points 3, 7, ami 1 , rcspcttively, in f ifjurc S. It is normal for the 
water of rivers entering lake Michigan and Cirecn Bay tu be warmer than the lake water during the 
spring and c,irly summer seasons, and dil ferentes as high as 20 T were measured by the investigators. 
Therefore, it is diflitult lif not impossible) to estimate luiw mucli of the relative warmth ol these 
rivers is due to natural kaiises and how muih is from industrial sources. Low-altitude, nighttime 
thermal 5t.tns would provide much more information on this type of problem. 

A noticeable aiiuatic thermal pattern is illustrated at point 2 (f ig. S). This pattern is interpreted 
tu be related to major circulation systems within lirccn Bay, as the confincmcnl of the pattern to 
the water area suggests that litc signal did not result from atmospheric tcmpcr.iiurc differences, ihc 
original RS>7 imagery provided numerous pallcrns related to lake currents and long shore currents 
in lake Michigan. The elfeci was similar to current turbidity patterns seen in die s.imc areas on 
IRTS MSS blinds 4 and 5 at approximately the same scale as the RS-7 imagery. 

Therm, vl patterns on tlie land surface were complex. The most readily dislinguishabic patterns 
were found over forested areas (for example, point 7, f igurc 5) which contrasted with those from 
cleared .igricultural areas and small lakes. Urban areas were nut distinct, but always gave a light lone 
on the original imagery, suggesting a "heat island" clfcct. An example is the city of Sheboygan 
(point 5). Tonal differences of lire urban areas compared lo the surrounding countryside resemble 
ERT5 MSS band 7 in appearance. 

One of the most noticeable terrestrial patterns on the RS*7 imagery is a dark swath with a north* 
west to southcusl oricnUitinn, located between Green Bay and Lake Winnebago (point (3). The 
probiibly cause of this cold signal is judged to be precipitation occurring shortiy before the time of 
fligin, resulting in higher soil moisture content in litis area. The only weather stations located within 
the dark'toned area are Brillion and Kewaunee, which reported 0.68 and 0.69 in. of precipitation, 
respectively, on 4 June (see App. A). 

RS-/4 thermal imagery. R5-1 4 imagery was obtained during Missions 272 and 273, which were 
flown on 31 M.iy and 12 June 1974 at approximately 2100 to 2400 hours local time (flighl lines 
are indicated in Figure 3). RS-14 imagery was also attempted for Mission 235 (June 1972) but 
ground fog obscured the results. Aircraft coverage of the Fox River was at 2500 ft (radiometric 
altitude); the rest of the coverage, including tliat of the Lake Micliigan shoreline and the Manitowoc, 
Sheboygan and Milwaukee Rivers was at approximately 3500 ft (radiometric altitude). This study 
focused primarily on the RS-1 4 imagery of the Manitowoc River and briefly on die larger scale cov* 
crago of the Fox River, Daui resulting from the PRT-5 radiometric thermometer intended for cor* 
relation with the RS-1 4 imagery were not received from NASA. 

The RS-14 thermal scan provided imagery in the 8* lo 14*pm wavelengths. The imagery was 
acquired at night so that only emitted thermal radiation would be recorded. The clarity of the 
thermal imagery acquired during Mission 272 is excellent. Except for the presence of scan lines, 
much of the imagery approaches photographic quality in terms of information content available for 
interpretation. 

The primary consideration in this study was the application of RS-14 thermal imagery to the 
Corps of Engineers mission work. The acquisition of thermal imagery is considerably more expen- 
sive than that of photographic imagery. This discussion compares the information content of the 
thermal imagery with color infrared photographic imagery acquired during Mission 273. 

Details of the RS-14 thermal scan of the Manitowoc River were compared with a 1:60,000 scale 
color infrared photograph acquired 23 days prior to the ilicrmal scan (7 May). Figure 6 is a nega- 
tive rendition of this RS-14 Imagery. Very little water detail Is visible In either type of imagery 
except at the wider portions near the river mouth. On the RS-1 4 imagery, temperature differences 
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on ihc RS*I<1 iniaRcry did nol create an idcniiriahlc signature on either higli altitude tnlnr nr color 
infrared imagery. In the study area> color infrared and/or color imagery acquired at a similar scale 
to the thermal Imagery could be used to locate most types of cflluents. I lowcver, details of ^lti< 
plumes, such as the si^e and lonccntratiun within mixing zones, w«juld be far better depicted by 
thermal imagery. 

Conclusions. Althougit much detail on land use patterns is portrayed in the tlicrmal imagery, 
its application In a mapping program siiould be limited to a few specialized areas, f or general map* 
ping cfforls the low icsolutiun, narrow Held of coverage, lack of geometric control and iiigh cy,< 
pense preclude the use ot thermal imagery. Por speciili/cd applications, however, the thermal scan- 
ner may be the only feasible sensor; for example, the dclination of impervious areas within a city, 
or differentiation of small drainage patterns witiiin a wetlands area must be done with a scanner. 

PM IS imagery 

On 31 May and 24 June 1973, the NASA NP-3 aircraft flew the RS-14 infrared scanner and the 
passive microwave imaging system (PMIS) over tlic test area. The PMIS is a recently developed sen- 
sor which has not been tested for many potential remote sensing applications. A variety of micro- 
wave sensors have been used to measure parameters of water or ice surlaccs, and several federal 
agencies have sponsored successful rcscarcli on determination of sea state, water surface temperature, 
water salinity, oil pollution of water surfaces, and sea mapping. 

PMIS system characteristics. Passive microwave systems do not provide tlicir own energy source. 
Tlicse sensors respond to the energy emitted and reflected from surfaces scanned by the sensor. 

The sum of tlicse types of energy constitutes the hrighlncss temperature of an object, fimitted and 
rellcctcd energy have an inverse relationship, with their sum being equal to one. Thus, objects which 
arc good reflectors of energy, such as metal roofs and areas of concrelo, are generally poor emitters 
of energy and vice versa. 

The greater variance of reflectivity in any given area, as opposed to variance In cmissivity, causes 
reflected energy to make up the major portion of a microwave signal recorded during daylight hours. 
Nighttime imagery records mostly emitted energy. This variance in reflection and emission causes 
the differences perceivable on the imagery. According to the Cartli Observations Aircraft Remote 
Sensing Handbook prepared by NASA (1972), the airborne PMIS gathers brightness temperature 
data (10,69 GHz, in keivins) from two scanning antennae, one for vertical polarization and one for 
horizontal polarization. Black and white images or. a video monitor and digital magnetic tapes arc 
generated aboard tlie aircraft. 

The imagery supplied for evaluation in this study was a computer-processed false color rendition 
of the scene scanned by the PMIS, each color representing a temperature range of about 4 K. A 
total of 64 colors were used, A computer printout of the area represented by the color transparency 
was supplied with the temperature for each footprint (color segment) given. 

The PMIS imagery was acquired on a single flight line of the NASA NP-3 aii craft at an average 
altitude above the terrain of 4300 ft. The flight line was centered on the Fox River, which drains 
Lake Winnebago and discharges into Green Bay about 40 miles to the northeast. The land use 
patterns along the Fox River covered by the PMIS swath include most of the land use categories to 
be found in the eastern Wisconsin test area and form a good basis for examining the applicability 
of PMIS imagery to patterns in the Great Lakes region. Recent studios conducted at the University 
of Cilifornia, Santa Barbara, have suggested dial passive microwave imagery' could bo useful for 
general land use mapping. 

The PMIS system was found in this study, however, to luvc severe limitations as a remote sensing 
tool for land use mapping. These limitations include low resolution and small area coverage. The 
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P~ location of aircraft (scanner) 
c = ground position, or center of scanner's field of 
viev/ 

a = altitude of aircraft above that of footprint 
N ~ nadir of aircraft at elevation of footprint 
G - horizontal distance from N to near edge of foot* 
print along Nc 
= a tan (ip~OI2) 

= 1.1180 

C = horizontal lengtii of footprint in direction of Nc 

= 0.108(7 
0=49.5“ 

ojj = vertical angular field of view = 2.6“ 

= horizontal angular field of view = 1.6“ 
^=£7sec49.5“ = 1.539(7 

w-2^ tan~j =0.0430(7 

Figure 9, PMIS scanning geometry, and formulas for footprint determination. 

scanning geometry of the PMIS system is shown in Figure 8, and the formula for determining the 
size of the PMIS footprint is shown in Figure 9. Applying the formula to the imagery evaluated in 
this report, the following dimensions are obtained: footprint length, 464 ft; footprint width, 185 ft; 
aircraft nadir to footprint, 5035 ft. It is interesting to note that from an altitude of 4300 ft the 
PMIS footprint size is slightly narrower but considerably longer than the 264-ft circle which com- 
poses a single unit (pixel) for the ERTS multispectral scanner from an altitude of 496 miles. The 
PMIS imagery evaluated has somewhat less resolution than the satellite MSS sensor. 

Several Interpretative approaches were used in the evaluation of the PMIS imagery. The initial 
approach was to visually match the major patterns of ^hc PMIS imagery with those seen on Aero- 
chrome IR transparencies and black and white prints of the study area. A Bausch and Lomb Zoom 
Transfer Scope (ZTS) was used to match the scales of the Imagery for visual comparison and inter- 
pretation. It was anticipated that major patterns identified on the photography as forest or 
agricultural land, or as residential oi industrial areas, would correlate with the microwave imagery. 
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Table If. Typical I’MIS signatures. 


Impmturc 


rantje 

JbL.. 

50-82 

Polarization 

Interpreted plwnomena 

Vertical 

Horizuntal 


83-122 

Vertical 

Hori/ontal 

Water, agricultural land 

123-226 

Vertical 

Hori/ontal 

Water, wetlands, industrial areas (metal roofs) 

Water, wetlands, residential areas, industrial areas (metal rnols) 

227-258 

Vertical 

Hori/ontal 

Boat berths, residential are.rs, industrial areas (metal riinls), 
agricultural areas, heavily vcgctateil areas, isolated mcial- 
runfed buildings, trailer parks, bare soil areas 
Rosidcnlial areas, Industrial areas (metal roofs), agricultural 
areas, heavily vegetated areas 

259-300 

Vertical 
1 lori/oiual 

Industrial areas, residential areas 
Residential areas 


It was not possible to identify the major patterns on the PM IS imagery. The only dependable 
p.ntcrn agreement was with the land-water boundaries representing the Fox River and the shore of 
Green Bay. In some cases other kinds of differentiation could be observed, for instance, the boundary 
between the central business district in the city of Green Bay and the established residential areas with 
many trees. Wet areas and small ponds produced a distinctive signature on the PMIS imagery and 
could ho compared to the photo record for shape and si/c. However, such relationships were not 
consistent from place to place, and pattern interpretation based on the PMIS Imagery without con- 
ventional photography as ground truth was not feasible. 

An attempt was made to correlate the PMIS computer printouts to color transparencies, A direct 
correlation could not be made, and in many cases, a unique footprint shown in the color transparency 
could not be located on the digital printout for the corresponding scan line. Nevertheless, it was 
possible in most cases to correlate a particular cultural or natural feature as shown on the color infra- 
red with outstanding footprints on the PMIS imagery. A summary of interpreted relationships is 
sliown in Table II. Possible causes of problems in correlation are: 1 ) individual scan lines are not 
properly syndirtini^ed with the computer printout, 2) the imagery may show unique segments when 
the brightness lemporaiures do not change significantly, or 3) the lime index printed on the imagery 
or the lime indication on the computer printout may be incorrect. The best method of evaluation 
was a simple visual comparison of the PMIS Imagery to conventional photography. 

Cvoluadon of PMIS imagery. The particular PMIS computer imagery cv.iluaicd in this project 
did not provide satisfactory results for the purposes of investigating land use. Wetlands wore visible 
on the imagery, but the possibility of confusing such terrain with various other phenomena that pro- 
duce similar tonal signatures was too great to be overlooked. One cluster of buildings would appear 
as a white or pink area on the imagery, while a similar group of adjacent buildings would n'j)l be repre- 
scnled at all. Of two adjacent ponds of similar size, one would be visible and the other not shown. 

On the basis of color or tone signature, heavily vogcialed areas could not be scp,iralccl from residen- 
tial or industrial areas. In some cases, land could not be distinguished from water on the microwave 
Imagery. Some roofs showing white on the imagery could be mistaken for ponds, as water also may 
be white on the imagery. This type of identification problem dominated corr<'lativc attempts based 
upon color and lone on the PMIS imagery. 
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a. t>. 

f igure tO. Comparison of high-attitude color infrared imagery (a) 
with PXilS imagery (b). 
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Wilh regard lu possible practical applicalions of the PMiS as a remote sensing tool in the eastern 
Wisconsin test area, it virould appear there arc severe limitations, especially for terrestrial applications. 
These limitations include low resolution, smalt area coverage, difficulties of scan line identification, 
and problems with numerical correlation of the colors represented on the imagery. This negative 
evaluation did not, however^ apply to the traditional applications of passive microwave sensors which 
arc primarily related to water surfaces.* The few areas of wetlands or fill were lor the most part 
detectable, A comparison of PMIS imagery with RC-S photographic Imagery Is shown in figure 10, 
The scale of the PMIS imagery in Figure 10 is slightly smaller in its long axis compared to the RC-8 
Photographic image, so that comparable points do not register exactly left and right. Although both 
originals were in color, this bla k and white reproduction provides considerable information on sig- 
nature comparability. The major water bodies, Green Bay and the Fox River, are easily defined in 
both types of imagery. The residential areas (point 4, Fig, 1 0) appear as a warm (dark) signal in liic 
PMIS, probably due to dense tree growth. The central business district of Green Bay (point 2) gives 
a mixed signal. The arc of cold (light) signals just below this point is an error, due to a tilt of the 
sensing aircraft at this point. Paper mills (point 3) provide a cold PMIS signal, and arc highly reflec- 
tive in the photographic wavelengths, 

f or land use mapping and other terrestrial applications in the eastern Wisconsin test area, the 
PMIS docs not appear to be a practical tool. The only clear advantage it has over other available sen- 
sors is that imagery can be acquired under adverse weather conditions. Although more delailcd pro- 
cessing at tlhe PMIS ground data station and imagery acquired at a lower altitude would provide more 
detail, it appears unlikely that the PMIS could provide data for land use mapping comparable, for 
most practical applications, to cither thermal scanning or photographic sensors. 


LAND USE ANALYSIS 

The initial step in the land use analysis portion of this study was to determine mapping criteria 
which would permit a quantitative comparison to water quality measurements. Most urban and 
many agricultural point sources can be directly related to a potential for pollution, but meaningful 
areal measurements cannot be obtained. Although changes in number and types of point sources 
will occur as land use patterns evolve, the size of ihe site or land area occupied by a typical point 
source is not directly related to the intensity or type of pollution. The developmentof a data base 
for land usc/watcr quality comparisons was, therefore, restricted to land use categories that lend 
themselves to area measurements rather than enumeration of point sources. This investigation was 
then focused on the regional relationship of land use categories to sediment loading in streams. 

The uses Land Use Classification System was selected for use in this investigation. The benefits 
of this system arc described (Anderson et aL 1972) as: 

The classification system utilizes the best features of existing widely used classifica- 
tion systems to the extent that tlicy are amenable to use with remote-sensing, and it is 
open-ended so that regional, state, and local agencies may develop more delailcd land 
use classification systems, at third and fourth levels, to meet their particular needs and 
at the same time remain compatible with each other and with the national system. 

An outline of the system is shown in Table IIL This system is adaptable to three levels of organiza- 
tion. Levels I and II arc readily mapped from most types of aircraft photography and require a minf 
mum of ground truth data for verification. The only extensive ground truth needed would be that 
required for some level 111 land use mapping. Level II mapping was accomplished on a small area 


*The use of PMIS for the differentiation of soil moisture was not sufficiently tested, due to the 
uniformity of soil type within the Fox River valley, 
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primarily to compare conient and efficiency of mapping from tRTS-1 vs high allituclc aircraft 
imagery (Haugen and Spicll 1974). The only level 111* categories utilised were to diffcreniiale be- 
tween open agricultural land and dial with a cover crop. 

For land use mappings the color infrared imagery was primarily acquired from the RC-S camera 
(6-in. lens), with an approximate scale of 1:1 20,000. (Although Zeiss 1 :60,000 photography pro- 
vided the greatest detail, mission rcslriclions precluded complete area coverage with tills camera.) 

The RC-8 color Infrared transparencies were contact printed to obtain black and while negatives and 
were enlarged to a print scale of 1 :63,360 for use with acetate watershed outline maps prepared from 
uses 15-minute topographic maps. Inlcrprcialion was done directly on the acetate overlays. These 

9-in. color transparencies were placed on a nearby light box for constant reference during the 
mapping on the black and while enlargements* 

Analysis of 1972 imagery 

The East and West Twin River watersheds were selected for the initial mapping effort because the 
photography showed obvious color-density differences between the two streams, v/hich suggested 
apparent variations in sediment load. The object was to correlate the apparent changes in sediment 
load to land use. Since the 1 ; 120,000 color infrared imagery provided the greatest amount of infor- 
mation on land use differences Jt was used for land use interpretation. Basic details such as roads, 
streams and watershed boundaries were transferred directly from the base map to an acetate overlay 
with predrawn linear features to provide control for distortions on the black and while photo en- 
largements. Colors were chosen for the manuscript map, based on their optical density, to allow 
measurement of areas through color densitometer planimetry. The photo base map and the manu- 
script land use map arc illustrated in Figures 11 and 12. Unfortunately the colors most suitable for 
color densitometry do not reproduce well in black and white photography, so that the categories 
arc difficult to discern in the illustration. 

The inslrumcnl used was an Anlcch model A-1 2 densitometer, which consists of a black and white 
TV camera especially modified for uniform sensitivity across its entire scan area, a color TV receiver, 
and an electronics package including a digital planimctcr* The instrument converts density levels 
into different voltages which are then displayed as discrete colors and/or numerical values. The 
range of density values examined is adfuslablc, and the assignment of specific colors to given density 
values is operator-controlled. A "color window" feature permits the analysis of a variable-size square 
or rectangular area within the scanned area. The planimctcr feature provides a direct readout of the 
proportion of the window in a given color or density level* Using the densitometer planimetry 
method, 97-99% accuracy is achieved, provided the initial mapping is correct. 

The color window must be calibrated to correspond to a known area on the acclate or base map. 
Since the mapping was done at a scale of 1 In. = 1 mile, a grid network (Fig. 13) representing 16 in.^ 
on the base map, or about 10,240 acres on the ground, was established. This block, size was selected 
because it yielded maximum results in terms of accuracy and size measurement. The grid was ap- 
plied to the Twin Rivers watersheds and the resulting measurements arc indicated in Tables IV and V. 

In the above analysis an experiment was accomplished by using the color densitometer to measure 
land use areas directly from the color and color infrared imagery. Color filters were employed to 
emphasize tonal patterns and to vary contrast differences indigenous to the imagery. The obvious 
advantage of this procedure is that manual mapping of land use is eliminated and data for compari- 
son to stream sedimentation and other water quality measurements can be obtained directly. Also, 
maps of land use can be compiled by simply photographing the individuai sections on the densito- 
meter screen. 


^ The level 111 classification in this system is used to further subdivide level M categories. The choice 
of categories within level 111 is made by the investigator to fit the objectives of his study. 
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hqurv II Photo hasr nujp with watershed outlines of f ast and West Twin Rivers. 
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Figure 12. Bfaclf and white reproduction of manuxript land use map. 


The direct measurement of land use areas from the photography proved to be feasible only tor 
limited areas. Only thos^ ai^as with light-toned soils provided sufficient color density contrast be- 
tween open fields and vegetation. Other land use categories, notably forests and roads, almost al- 
ways could be distinguished and directly measured. However, the method was not considered feas- 
ible because all level I categories could noi be measured exactly from place to place within water- 
shed areas. 

The next approach to land use mapping involved hand-colored overlays on frosted acetate, but 
with three important differences from the original effort on the Twin Rivers watersheds. First, a 
subcategory was created to distinguish **open** agricultural land, essentially areas of bare soil where 
crops had not yet reached sufficient density (maturity) to produce a signature on the color infrared 
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Figure 13. Grid used far dens/toineter measurement. 


imagery. (All other agricultural land was designated as ‘‘closed.”) Open land, being more subject to 
erosion, would contribute more directly to turbidity differences in streams. The turbidity differ- 
ences were apparent in the imagery and could be measured by the analysis of water samples. Secondly, 
this mapping effort was done directly on an overlay of the 1 :1 20,000 color infrared imagery. The geo- 
metric control originally provided by USGS base maps was not believed critical to the measurement 
of land use areas, The smaller scale, together with the procedure of mapping directly from the color 
Infrared imagery, was considerab'; faster. The third difference was the production of a strip map 
which furnished a sampling of land use units representative of the watershed area. Two grids - one 
with two square miles, and one with eight square miles - were prepared for equal sampling to Include 
at least one mile on each side of the stream. The small grid square was used where the stream was 
more or less straight; the larger grid was used where streams bifurcated or where major curves were 
encountered. 
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Table lit. Land use classification system for use with 
remote sensor data. 


l evel I 


Lewi It 


01, Urban and built-up land 


02. Agricultural land 


03. Rangeland 


04. f-oreslland 


05, Water 


06. Nonforested wetland 

07. Barren land 


08. Tundra 

09. Permanent snow and teefieids 


01. Residential 

02. Commercial and services 

03. Industrial 
O'!. Lxtractivo 

05, Transportation, communications and utilities 

06, Institutional 

07, Strip and clustered settlement 

08, Mixed 

09, Open and other 

01. Cropland and pasture 

02. Orchards, groves, bush fruits, vineyards, and 
horticultural areas 

03. Feeding operations 

04. Other 

01. Grass 

02. Savannas (palmetto prairies) 

03. Chaparral 

04. Pcsert shrub 

01. Deciduous 

02. Evergreen (coniferous and other) 

03. Mixed 

01. Streams and waterways 

02. Lakes 

03. Reservoirs 

04. Bays and estuaries 

05. Other 

01, Vegetated 

02. Bare 

01. Sait flats 

02. Beaches 

03. Sand other than beaches 

04. Bare exposed rock 

05. Other 

01. Tundra 

01. Permanent snowficlds and icefields 


Two watersheds were mapped and measured using this procedure, the Oconto (1060 square miles) 
and the Manitowoc (557 square miles). The Oconto, located in the northern part of the test area, is 
in a predominantly forested area with few apparent sources of aquatic pollution. The Manitowoc 
watershed had a proportionment of land use categories typical of the eastern Wisconsin test area. 
Measurements of land use areas for the Oconto and Manitowoc watersheds were determined by use of 
the strip mapping method. This technique resulted in a mapping area approximately one-third of the 
total watershed. 
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Table IV. Area of land use types, East Twin River watershed 
{84,697 total acres). 


Measumd 
block 
(Fig, 13) 

Forest 

Urban 

Agriculture 

Other* 

9( Acres 

% Acres 

% 

Acres 

9i Acres 

1 

471 

1301 


860 

72 

10 

1608 

184 


4413 

420 

9 

850 

522 


2775 

102 

13 

503 

0 


1577 

123 

12 

2826 

0 


6093 

317 

11 

993 

0 


3052 

358 

20 

922 

358 


4547 

287 

19 

3686 

10 


5642 

481 

18 

707 

0 


2028 

225 

23 

236 

0 


358 

10 

22 

2458 

0 


5673 

113 

21 

1782 

0 


5243 

277 

30 

1341 

0 


4024 

440 

29 

3389 

0 


5335 

543 

28 

348 

0 


778 

51 

32 

1065 

0 


635 

307 

31 

266 

0 


297 

113 

Totals 

28 23,452 

3 2375 

64 

53,330 

5 5,257 84.697 


* This category includes watershed boundary lines, roads, rivers, etc. 


There is a bias in the use of strip mapping for measurement of land use categories within a water- 
shed, because land use within two miles of a stream is not necessarily representative of the propor- 
tions of land use within the entire watershed. However, for purposes of studying the relationships 
among land use practices, stream sodimontalion and other water quality measurements, it can be 
argued that the resultant data from this procedure are pertinent. 

Analysis of 1974 imagery 

Several changes in the mapping approach were made, based on experience gained in utilizing the 
1972 imagery for the mapping effort described above, instead of using colors on a single overlay 
map as in the initial procedure, greater measurement accuracy could be achieved with the densito- 
meter/planimeter by delineating land use with black ink on separate overlays For each category. Al- 
though mapping with colors is the more rapid method, boundaries are not as distinct and measure- 
ment with the densitometer is more difficult. This is due to the unavoidable density differences in 
the colored patterns on acetate and also to the similarity of color densities. 

The East Twin River watershed was selected as a test area for this mapping approach. This water- 
shed is small (approximately 131 square miles) compared to the Oconto and Manitowoc watersheds, 
and water samples were obtained at several locations along the East Twin River. The 1974 color 
imagery was used to prepare a black and white photographic base map at a scale of 1 :63,360. The 
1974 color infrared imagery was overexposed and could not be utilized for enlargement. 

Five acetate overlay maps were prepared for the East Twin River watershed at a 1 :63,360 scale 
representing agriculture (open), agriculture (closed), forest, urban, and water. This total watershed 
map provided an accurate basis for evaluation of previous and future mapping techniques. A repro- 
duction of the base photo map together with overlay for agriculture (“open") is shown in Figure 14. 
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Table V. Area of land use types, West Twin River watershed 
(106,205 total acres). 


Measured 

block 

Forest 

Urban 

Agriculture 

Other* 

96 

Acres 

% Acres 

% Acres _ % 

Acres 

1 


1608 

1638 

1556 

625 

2 


1925 

430 

4076 

358 

3 


922 

0 

2252 

717 

7 


696 

287 

2396 

819 

15 


819 

0 

2488 

297 

14 


1137 

573 

6154 

225 

S 


1751 

522 

7393 

522 

9 


655 

10 

3553 

388 

12 


614 

0 

102 

10 

13 


2990 

92 

4014 

625 

16 


1864 

0 

7905 

205 

17 


2171 

840 

6728 

317 

18 


1475 

0 

4618 

164 

25 


666 

0 

3553 

143 

24 


1577 

0 

6738 

1034 

23 


2488 

0 

5612 

522 

22 


20 

0 

789 

0 

28 


963 

0 

1864 

92 

29 


41 

0 

133 

0 

27 


819 

369 

3717 

72 

26 


584 

0 

789 

225 



0 

0 

92 

0 

33 


215 

0 

1382 

0 

Totals 

22 

16,000 

4 4761 

67 77 904 6 

7,360 106,025 


* This category includes watershed boundary linos, roads, rives, etc. 


The time required to map an entire watershed at the 1 :63,360 scale was too great with regard to the 
principal obicettve of obtaining quantitative land use data for several representative watersheds. 

Therefore, a second sampling technique was devised and tested for precision against the totally 
mapped East Twin River watershed. One square mile (one section) was established as the sampling 
unit. A computer-generated listing of random numbers was obtained to provide a 25% sample for 
grid matrices of 3x3, 4x 4, 6x6, and 12x 12 units. The different sizes were necessary to provide a 
"best fit" when the grids were systematically overlaid on the irregularly shaped watershed. The grid 
is illustrated in Figure 1 5, The grids were applied to a watershed mosaic composed of USGS topo- 
graphic maps. The random numbers, identifiod as to row and column for each grid, were used to 
locate the square-mile sample areas, Frequently roads, conspicuous on the Imagery, followed section 
lines. The exact location of a section-sized area was varied slightly from the location indicated by 
the grid, so as to conform to these boundaries for ease in subsequent mapping. For the East Twin 
River watershed the 25% sample provided 36 one-mile squares to be mapped. 

Measurement was done with the color densitometer as previously described. In this case, two 
factors contributed to greater measurement accuracy; 1 ) the color window could be calibrated to 
conform exactly to the one-square-milc size, and 2) the black Ink patterns created only one density 
level so that the density settings of the instrument remained constant. To verify the sampling method, 
all East Twin River watershed map overlays were measured using larger grid squares and the results 
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/ iqurr 14, Photo mouiu of i os/ Iwm River wurenhed with {a) hose map and (b) ovtrhv for 

agncutturat (dosed) fund. 


were ciimp^red. Ihcsc results indicated that the 25% sample was within 8% of the values obtained 
from the total sample for all categories and within 4% for the maiority of categories. The bias ex* 
pericneed with the strip mapping approach was avoided. 

f or the mapping and measurement of the OconU) and Manitowoc watersheds, the 25% sample 
approach was also used. Overlay maps were produced from the 1 : 1 20,000 color infrared >mager 
This imagery, although over-exposed, contained enough detail to permit direct iiiterpretation. 
f unher timesaving resulted from combining the mapped categories on a single sheet of acetate. 
Seven separate borderless squares, representing each of the categories, were mapped in columnar 
ft>rm on an acetate overlay. This permitted mapping alt factors within the section sequentially, and 
provided a systematic basis for the subsequent densitometer measurement of areas. In this manner, 
level I mapping of a square-mile section could be accomplished in approximately 7 to 1 5 min, de- 
pending or the number of categories present and pattern complexity. An illustration of these over- 
lays is sh( wn in f igyre 16. 

SOILS OF MANITOWOC. EAST TWIN AND OCONTO RIVER WATERSHEDS 


Assessment of the soils was needed to examine the relationship among land use, soil association 
and sueam sedimentation. The soil information used in this investigation was obtained from soil 
surveys and maps prepared by the U.S. Department of Agriculture. Deuils on soil association, 
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Ongmol scate 1m ^ Irm 

Figure J 5. Grid matrices for a 25% random sample of variously sired 
mapping units. 


Oconto River 
Sample No. 204 


Figure 16. Example of separate overlay mopping for a single square- 
mile section. Numbers refer to land use categories. 



topograpiiy, parunt material and classincation (according to the 7lh approximation) can be obtained 
from the legend of the Overlay Soil Map of Wisconsin (Wisconsin Geological and Natural History Sur* 
vcy 1968). 

Manitowoc River watershed 

The approximate percentages of tlie principal soil groups that occur In the Manitowoc River 
watershed (Fig, 17} arc: soils of the Southeastern Upland (20%), soils of the Northern and Eastern 
Clayey and Loamy Reddish Drift Uplands and Plains (50%) and soils of the Stream Bottoms and Major 
Wetlands (30%). The Southeastern Upland soils, developed in loess, sand or loamy till, arc located 
in the headwater area of the watershed. The parent material for the soils found on the Reddish Drift 
Uplands and Plains is primarily clay and silt. These soils occur in the central and lower reaches of 
the drainage system on 2% to 12% slopes. The Stream Bottom and Wetland Soils arc located through- 
out the watershed on level plains and arc usually poorly drained. The parent materials for these soils 
are glacial sandy oulwash and organic materials, calcareous clayey lacustrine sediments and sandy and 
loamy glacial drift. 
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Figure 1 7. Soil groups of the Marti towoc River watershed. 


East Twin River walershcti 

The same soils groups that appear in the Manitowoc watershed occur in the East Twin River 
Watershed {Fig. 18), but the percentages differ asTollows; soils of the Southeastern Upland (55%), 
soils of the Northern and Eastern Clayey and Loamy Reddish Drift Uplands and Plains (40%) and 
soils of the Stream Bottoms and Majoi Wetlands (5%). The distribution of these soils is similar to that 
of the Manitowoc. 


Oconto River watershed 

The Oconto Is the largest of the three watersheds and contains the greatest number of soils. The 
appro^dmatc percentages of the principal soil groups (Fig, 19) arc: soils of the Southeastern Upland 
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Ffgi/re 18 . Soil groups of the East Twin River watershed. 


and Plains (5%), Solis of the Northern Loamy Upland and Plains (22%), soils of Northern Sandy Up- 
land and Plains (1 7%), and soils of Stream Bottoms and Major Wetlands (7%). Generally the silty and 
loamy soils occur on moderate to steep slopes in the headwater regions. Coarse-grained deposits are 
found in the gently rolling and level plain positions in the central portion of the watershed. Finer 
textured soils (silts and clays) arc found in the lower reaches of the river. 
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Figure 19, Sol! groups of the Oconto River watershed. 
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WATER QUALITY DATA 


Although the primary source of information on water quality throughout titc study was laboratory 
analysis of water samples collected in the field, efforts were also made to derive water quality infor- 
mation by remote sensing means. The initial attempt was the attempt to develop a water color classi- 
fication based on stream colors observable in the 1972 1 ;1 20,000 scale Ektachrome (SO-397) imagery. 
Some of the shortcomings of this imagery for the purposes of pollution Identification have already 
been discussed. Its small scale, the effects of atmospheric attenuation, and the loss of contrast and 
detail due to two generations of reproduction (internegative and duplicate positive transparency) of 
the imagery supplied by NASA all serve to lessen the amount of information which may be derived. 

It was decided, however, that sufficient color detail and uniformity existed in this imagery to make 
possible comparisons for land use classification. 

Water mass color determinations were done according to the Densitomcter-Munsell system of 
color notation (Rib 1968). This technique permits the determination of the Munscll notation with 
a transmittance densitometer so that errors involved with matching of color chips are eliminated. 

The measurements were made with a Macbeth TD-504 Quantalog transmission densitometer, fitted 
with a custom-made 0.75-mm aperture so that measurements could be made of some of the smaller 
streams. The instrument is normally supplied with 1- and 2-mm apertures. The measurement tech- 
nique for each point consists of separate readings through each of four Wratten filters: 106W (visual), 
92 (red), 93 (green), and 94 (blue). The raw measurements are converted to Munsell values of hue, 
chroma, and intensity according to the procedures outlined by Rib (1968). 

The feasibility of correlating tonal variations of aerial imagery w[»h turbidity measurements for 
local areas has been demonstrated (Shorz et al. 1969 and Lillesand 1975), The effort In this study 
was to attempt a more general relationship that could be applied, on a regional basis, primarily for 
correlation with agricultural land use. The results of numerous measurements of color signatures of 
streams throughout the study area indicated that this approach was probably not feasible. Although 
there was sufficient uniformity of color representation of the maior streams within the 1972 color 
imagery, the color balance of the 1 974 imagery of the same locations was completely different. In 
the comparison of streams where gross differences in the level of turbidity or suspended solids were 
known to exist (for example, the Fox as compared to the Oconto Rivers), differences could be seen 
on all the imagery. But in general, it had to be concluded that the representation of stream color In 
high altitude aerial photography displayed far more variation due to factors other than turbidity in 
the streams than it did f.n' these parameters. 

To provide a more direct basis for comparing water quality parameters to signatures contained in 
remote sensing imagery, multispcctral measurements of refieclcd and incident light were made in 
conjunction with water sampling during Mission 235. A non-imaging hand-held four-band Exotech 
ERTS Ground Truth Radiometer was employed. 

Multiple correlation analysis indicated that reflectivity of the streams sampled is related linearly 
to incident and reflected radiation In ERTS MSS band 5 (Fig, 20), A multiple correlation coefficient 
R of 0.85 and a standard error of estimate for suspended solids in streams of 3,8 mg/1 was obtained. 
The measurement and estimated values are shown in Table VI. The range of suspended solid concen- 
trations during the test period was low, 2*20 mg/I, but the high correlation with radiated and incident 
light measurements suggests that the relationship would also exist for higher concenfations of sus- 
pended sediment in the streams. Further sampling of a wider range of suspended sediment concen- 
trations with concurrent incident and radiated light measurements would permit the estimation of 
stream sediment loads based on digital density measurements of remote sensing imagery In the 0.6- to 
0.7-/um wavelength for this test area. 
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Figure 20. Relationship of E RTS MSS band 5 Incident and reflected 
radiation to measured suspended solids. 

Table VI. Comparison of spectral measurements* and suspended 
sediment load at water sampling sites. 


Suspended solids 

Stream tSampHng Date! Reflected radiation band Incident radiation bandf (mglD 

me Local time J 2 3 4 1 2 3 4 Meas. Est.** 


Oconto River 

2 May 1974 

0.04 

0.05 

0.03 

0.00 

0.30 

0.38 

0.30 

0.34 

3.4 

2.2 

Hwy 141 

1530 











Oconto River 

2 May 1974 

0.10 

0.12 

0.08 

0.04 

0.30 

0.40 

0.32 

0.36 

1.4 

6.2 

Stiles 

1600 











Manitowoc River 

3 May 1974 

0.12 

0.20 

0.10 

0.02 

0.58 

0.70 

0.18 

0.08 

10.3 

6.9 

Clark's Mills 

1120 











Manitowoc River 

3 May 1974 

0.22 

0.32 

0.18 

0.08 

0.66 

0.82 

0.64 

0.64 

16.9 

12.6 

Mill Road 












Manitowoc River 

3 May 1974 

0.30 

0.42 

0.20 

0.08 

0,56 

0.80 

0.64 

0.80 

16.5 

19,2 

Harbor 

1300 











East Twin River 

3 May 1974 

0.14 

0.42 

0.14 

0.06 

0.60 

0.74 

0.60 

0.78 

4.1 

8.8 

Mishicot 

1330 











Fox River 

3 May 1974 

0.16 

0.24 

0.14 

0.00 

0.16 

0.16 

0.16 

0.60 

19.1 

17.2 

Wrightsiown 

1030 











Fox River 

2 May 1974 

0.16 

0.20 

0.11 

0.00 

0.28 

0.32 

0.24 

0.28 

12.2 

12.4 

Green Bay 

1315 











East Twin River 

3 May 1974 

0.16 

0.20 

0.08 

0.02 

0.46 

0.60 

0.46 

0.56 

10.1 

8.3 

Choto 

1500 











♦Spectral range 

of bands: 1 (0*5-0 j 

Sym); 

2(0.6 

•0.7 pm); 3 (0 

I.7-0.8 

pm);-! 

1 (0.8-1 .' 

I pm). 



fRcflcctance measurements are based on a 2*niW fuil scale reading (2 mW/cm^ ); incident measure- 
ments arc based on a 10-mW full scale reading (10 mW/cm^}. 

♦♦Estimated according to y - 4.54^62.5X-^4.5X^, whercy = estimated suspended solids (mg/1), A', = 
reflected 0.6-0.7 ym radiation (mW/cm*x5 ), Xz = incident 0.6-0.7 ym radiation (mW/cm*). 
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Water samples were taken in June and October 1973 and May 1974 as part of the ground truth 
for the NASA photo mission. Suspended solids were determined atCRREL according to procedures 
used by the Wisconsin Department of Natural Resources, Division of Environmental Protection and 
the U.S. Department of Agriculture, A location map of the field sample sites is shown In Figure 21 . 

Manitowoc River 

The June 1973 Manitowoc River suspended solids data varied from 30 mg/1 to 105 mg/i with an 
average of 70.1 mg/l (Table VII). The total precipitation, measured at local testing sites immediately 
preceding and during the sampling period, 19-21 June, averaged 0*45 in, (see App* A). During 
October 1973 the suspended solids in the Manitowoc River decreased to an average of 6.3 mg/l* 

The samples were obtained on 8 October, and for the six days, 3-8 October, only 0*1 to 0.2 in. of 
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Table VII. Average daily water quality data: Manitowoc River, 
June and October 1973 and May 1974. 


Sampling 

site 

Turbidity 

(ITU)* 

Suspended 
solids 
(mg II) 

Organic 

matter 

(mgll) 

Nonorganic 

{%) 

Organic 

(%) 

June 1973 







Clark Mills (2) 

43 


81.8 

25.4 

69 

31 

Mill Road (3) 

39 


29.8 

9.6 

69 

32 

Clark Mills (2) 

18 


63.4 

18.6 

71 

29 

Mill Road (3) 

52 


105.4 

16.2 

85 

15 



Avg 

70.1 




October 1973 







Hwy, 141 (3) 

2 


5.0 




Harbor (1) 

5 


7.6 






Avg 

6.3 




May 1974 







Harbor (1 ) 

10 


16,5 




Mill Road (3) 

10 


16.9 




Clark Mills (2) 

2 


10.3 






Avg 

14.3 





JTU - Jackson Turbidity Units. 


precipitation occurred in the watershed, in May of 1974 the average amount of suspended solids 
in the river was 14.3 mg/I. For the six days iu and including the sampling date of 3 May the 
precipitation averaged 0.2 in. The suspended iolids data for the same monthly time periods from 
1965-1968 indicate the same trend as the dati obtained from this study. 

East Twin River 

Suspended solids determinations were made at four sample sites in the East Twin River (Fig. 20). 
The average suspended solid concentration from the 21 June 1973 sampling was 35.4 mg/I (Table 
VIII) and precipitation during the preceding seven days was 0.39 in. The October data were collected 
on the 8lh when the suspended solids averaged 13.0 mg/I and precipitation for 3-7 October was 0.28 
in. The May sampling, taken on 3 May, yielded 4.1 mg/I of suspended solids. The precipitation 
totaled 0.22 in. for the six days from 28 April to 3 May. These data compared favorably with those 
obtained by the Wisconsin Department of Natural Resources (1965-1968). 

Oconto River 

The three water sampling sites for the Oconto River arc shown in Figure 20. The 19 June 1973 
samples yielded an average of 61.9 mg/I of suspended solids (Table IX). The total precipitation from 
all measuring sites during 14-18 June was 0.15 in. In October 1973 and May 1974 the water 
samples collected contained very little suspended solids and much of this was in organic form. 

These data correspond to data taken during similar intervals in 1965-1969 (Wisconsin Dept, of 
Natural Resources 1968). 
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Table VIII, Average water quality data: East Twin River, 
June and October 1973 and May 1974, 


Sampling 

site 

Turbidity 

(]TU)* 

Suspended 

solids 

(mil) 

Organic 

matter 

(mgfl) 

Nonorganic 

(%) 

Organic 

(%) 

June 1973 






MIshicoi (5) 

12 





Tisch Mills (6) 

18 

28,6 

8.6 

70 

30 

Two Rivers (4) 

15 

39.6 

8.2 

79 

21 

Co, Rd. B (7) 

26 

38,0 

16.4 

57 

43 

Strangoville 


Avg 35,4 




October 1973 






Hwy, 42 (4) 

20+ 

22.4 




Mishicot (5) 

4 

3.7 






Avg 13.0 




May 1974 






Mishicot (5) 

2 

4.08 




* Jackson Turbidity Units, 



Tfible IX. 

Average water quality data: Oconto River, 



June and October 1973 and May 1974. 





Suspended 

Organic 



Sampling 

Turbidity 

solids 

matter 

Nonorganic 

Organic 

sites 

(JTU) 

(mail) 

(mail) 

(%) 

(%) 

June 1973 






City 

24 

52.6 

9.2 

83 

17 

Hwy. 141 

35 

71,2 

16.8 

76 

24 



Avg 61.9 




October 1973 
Hwy, 141 
May 1974 

2 

2.2 

1.5 

59 

41 

Hwy. 141 

22 

3.4 




Stiles 

2 

1.4 






Avg 2.9 
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LAND USE/STREAM SEDIMENTATION RELATIONSHIPS 


The land use maps prepared from the NASA RB-57F imagery for the Manitowoc, East Twin 
River and Oconto watersheds Indicate that the agricultural land composes the majority of the land 
area in the Manitowoc and East Twin River, 74% and 66% respectively, but only 27% of the land 
in the Oconto watershed, (The malor unit in the Oconto watershed is forest, making up 50% of 
the landscape,) Witen these percentages arc compared to the suspended solids data taken from the 
Manitowoc, East Twin and Oconto River watersheds, the results show that as the percentage of 
agricultural land increases, the amount of sediment in the stream increases (Table X). 


Table X. Land use, suspcncied solid, soils and precipitation data for 
the Manitowoc, East Twin, and Oconto River watersheds. 


Watershed 

Total 

square 

miles 

Agr 

(open) 

(%) 

Agr 

(dosed) 

(%) 

Forest Urban 
(%) (%) 

Wetlands 

and 

water 

(%) 

Other 

(%) 

Suspended 

solids 

average 

(mti) 

Soil 

texture 

Total 
preclp 
five days 
preceding 
mission 
(in.) 

Manitowoc 

505 

35 

39 

17 

2 

5 

2 

14.6 

Fine 

0.2 

East Twin 

166 

36 

30 

IS 

5 

2 

9 

4.1 

Medium 
to fine 

0.2 

Oconto 

966 

18 

9 

50 

1 

16 

6 

2.4 

Coarse 

0.1 


Supporting data for this relationship also can be found in the soils data. The finer textured 
soils and the largest percentage of these soils can be found in the Manitowoc watershed. The head- 
water area of the East Twin contains medium- to fine-grained surface sediment, but toward the 
mouth of the river there arc finer grained soils. However, In the Oconto the soils arc primarily 
developed in sandy and loamy deposits which are very coarse-grained. Also, in this watershed 
there is a large percentage of organic soil, so that the suspended solids tend to contain a large 
percentage of organic matter. The Manitowoc River had the largest amount of particulate matter 
per liter during the test, These data reveal that the Manitowoc watershed has the greatest percentage 
of agricultural land and the finest textured soils. In conlr» »t, tb* Oconto has the smallest percentage 
of agricultural land, the coarsest textured soils and the lowest amount of suspended solids in the 
stream. 


RESULTS AND CONCLUSIONS 

The primary focus of this report has been to examine the feasibility of using remote sensing 
methods to rapidly and economically assess, on a regional scale, the effect of land use as it influences 
sediment loading of streams. A test area, consisting of several major v/atersheds in eastern Wisconsin, 
was selected for the development and evaluation of techniques to achieve this objective. 
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Sensor evaluation 

A varielv of aerial remote sensors were applied to the test area for evaluating and developing 
data pertinent Uj the assessment of Uie relationships between land use patterns and the sediment 
load in streams. By far the most useful imagery supplied by NASA for the investigation was the 
color Infrared photography acquired at 60,000 ft with 9-in. format RC-8 and Zeiss cameras. This 
combination of sensors was judged to be adequate for any regional scale land use mapping effort 
to level tl of the uses land use classification, and with appropriate ground truth support, level 111 
detail is feasible. 

Sensor evaluation included ilic following: 

High altitude RS*7 thermo! scanner. This imagery, acquired during daylight hours at an altitude 
of 60,000 ft provided a small scale (1; 1,000,000) thermal map of much of the study area. General 
land use patterns could be observed as well as large water circulation patterns in Lake Michigan and 
Green Bay. The detection and mapping of major circulation patterns in large water bodies would 
be the primary recommended application for this sensor for high altitude operation, 

RSH4 thermal scanner. This imagery, in the 8-14 pm range, was acquired for the V/isconsin 
shoreline and several watershed areas at altitudes range from 2500 to 3500 ft. It was compared for 
information content with color infrared photographic imagery obtained during other missions. 
Analysis indicated that many of the patterns displayed on the thermal imagery are readily detect- 
able on color infrared or color imagery. Recommended spcciatixcd applications, however, include 
monitoring and mapping details of effluent plumes for si^c and concentration, and differentiation 
within wetlands or impervious urban settings, or wherever relative temperature information is de- 
sired. For general mapping efforts, the lack of geometric control, narrow field of view, and high 
expense would appear to preclude the use of this sensor. 

High altitude multi-spectral 70-mm (Hasseiblad) photography. This sensing tool did not provide 
significant data input for the study. The small 70-mm format is difficult to work with, and the 
multispcciral and color infrared Imagery acquired with the Hassclblads did not contain any informa- 
tion that was not more easily obtained from the 9-in. format cameras. Although the acquisition 
cost of 70‘mm imagery is somewhat less than the cost of 9-in. format Imagery, inlcrpreiation for 
most purposes is more difficult and time-consuming and would probably offset any economy gained 
by the lower original cost. 

High attitude color (SO 397), This imagery was usable for most purposes of this study, but was 
surpassed in utility by the 9-in. format color infrared (2443) emulsion. Densitometric determinations 
of Munscll color equivalents were made at a large number of sample points. The results of this 
correlation effort were considered indecisive because of: 1) the wide variation in the general color 
balance of the imagery from flight to flight, 2) the loss of density detail on the duplicate imagery 
provided by NASA as compared to originals, and 3) the influence of many uncontrollable variables 
on photographic color, such as the angle and intensity of sunlight and atmospheric haze. It is 
believed that the method of correlating 50-397 imagery to water quality parameters could be 
developed to apply to small area, single flight line surveys, but it does not appear feasible for 
regional area surveys or with high altitude imagery. 

PMIS, The PM IS (Passive Microwave Imaging System) imagery was acquired on Mission 235 on 
the Fox River valley at an altitude of 4300 ft. Because PMIS is a new sensor, and few terrestrial 
applications have been demonstrated for it, considerable effort was devoted to the interpretation 
and evaluation of the imagery. Other than land-water boundaries, it was not possible to consistently 
identify major patterns on a 64 color photographic rendition of the PMIS scanning data. Further, 
exhaustive attempts to correlate Individual segments of the PMIS scan between the computer- 
generated false color imagery and the digital computer printout provided were unsuccessful, although 
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signatures of individual scan segments cuuld be correlated with photographic color infrared imagery 
for some types of features (such as large buildings with metal roofs or small impoundments in 
agricultural areas). Consequently there arc severe [imitations concerning the I’MIS as a remote 
sensing tool for terrestrial phenomena in a region such as the eastern Wisconsin test area. 

Land use analysis 

The principal objectives of the land use analysis portion of this investigation were to develop 
a data base reflecting regional differences within the tost area for comparison with stream sediment 
loading data, and to determine the most efficient and effective approaches to land use mapping 
based on the imagery provided by NASA. 

The uses Land Use Classification System for Remote Sensing was incorporated because it is 
widely used, is compatible with most other classifications, and is sufficiently flexible to accommo- 
date specialised applications. Although initial mapping efforts were confined to level I of the sys- 
tem, subsequent mapping incorporated levels It and 111 to provide more pertinent data on the 
erosion potential of agricultural land. Following the final mapping, comparisons were made with 
the earlier strip mapping efforts of the Oconto and Manitowoc watersheds and regional differences 
were noted. The strip map approach has a built-in bias which can vary, depending upon the land 
use development history and the nature of the terrain. The application of this mapping approach 
should bo confined to cases whore prior Information is obtainable on general distribution of land 
use types within the watershed. A flow diagram of the land use mapping effort is shown in 
Figure 22. 



Figure 22, Flow diagram of land use mapping and analysis. 
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Mapping approaches 

Various approaches to the actual mapping of land use were tested, but all incorporated the USG5 
classification system* The ability to measure mapped areas rapidly was a first requirement for any 
of the mapping procedures* The measurement device in all cases was the planimctcr component of 
a color densitometer. This required that all mapped categories be distinguishable on the basis of 
their gray scale density. 

The first large watershed map was done in color on frosted acetate and included only level 1 of 
the USGS classtfication, A concurrent attempt was also made to measure land use categories 
directly from color or color infrared imagery using the dcnsitomctric method. This method worked 
for level I mapping categories only *n limited areas with light toned soils. 

The second method applied to obtain land use data involved sampling rather than mapping the 
entire watershed. Colored acetate strip maps were made of two large watersheds, the Manitowoc 
and the Oconto, The sampling grid was constructed so that the sample included at least one mile 
of land on cither side of any stream. This made up approximately one-lhlrd of the entire watershed 
area. 

The third method involved a dual approach; one watershed, the East Twin River was measured 
first for the entire watershed area and then on the basis of a 25% random sample. The land use of 
this watershed was mapped in its entirety on six separate dear acetate overlays. The mapping 
was accomplished in Ink to simplify the densitomctric measurement and the results of this com* 
parison were that all the major mapping categories (Agriculture Closed, Agriculture Open, and 
Forest) were within 4 % of being equal. Larger differences were noted for categories such as Urban 
and Water where the area involved was quite small. Based on this comparison, it was decided that 
for subsequent mapping efforts the 25% random sample approach would provide sufficiently 
accurate data and result in greater efficiency. 

Land use/sedimentatlon relationships 

The analysis of the relationships between stream sedimentation and land use practices indicated 
tiiat the amount of suspended solids measured from water samples increased with the percentage of 
agricultural land within the watershed. Soils and precipitation data were also compared to the ^ 
suspended solids measured for the various streams. 

The use of high altitude color or color infrared imagery to characterize sediment loading of 
individual streams did not appear practical in the quantitative sense. The effects of Image duplica- 
ting, processing, photo color balance, atmospheric haze, and time of day of the imagery all served 
to decrease the informational content of the high altitude aerial images. Only major differences 
among the streams were observable, and these were evident only in some of the imagery. Attempts 
to measure and characterize stream color differences using the Densllometer-Munsell system of 
notation (Rib 1968) were not sufficiently consistent to be of value* However, a comparison of 
multispeclral measurements made with a hand-held radiometer in conjunction with water sampling 
during Mission 235 was successful, A multiple correlation coefficient/? of 0,85 was obtained be- 
tween incident and reflected radiation in the ERTS MSS band 5 {0.6-0.7 /im) and measured 
suspended solids concentrations in the streams. A trial application of the universal soil loss equa- 
tion (USDA 1965) did provide estimates of soil loss that correlated with observable and measurable 
differences in turbidity of the East and West Twin Rivers, but it was apparent that refinement was 
necessary for the input data to the universal soil loss equation to provide meaningful results for the 
variety of streams in the test area. 


40 


Cost effectiveness analysis 

A cost effectiveness analysis was perfornned on the various land use mapping approaches. The 
results arc summarized In Table XI. The strip mapping approach, using colored pencil on a single 
acetate overlay of original scale photography, provided the fastest method of obtaining land use 
measurements. Including analysis, 0.4 hours per square mile was required. For the areas sampled 
with a 2595 random sample and a combined sequential overlay map, 1.1 and 1.3 hours were spent 
per square mile of sample. The 25% sample of the East Twin River Watershed required slightly 
more time 3.2 hours per square mile of sample - because the mapping was done on separate 
overlays. The results of the cost effectiveness analysis are similar to those shown in a NASA- 
sponsored study by Vegas (1 973), although the mapping objectives were more general than this 
application. 


Table XI. Cost effectiveness analysis (in hours). 



East Twin 

Afanitowoc 

Oconto 

Area mapped 

25% 

Total 

25% 

25% 

Strip 

Enlargements 

1.5 

1.5 

2 

3 


Overlays 

4 

4 

8 

8 


Sampling 

Mapping 

8 


20 

20 

16 

Separate 

Combined 

29 

65 

120 

155 

40 

Measurement 

8 

12 

30 

35 

35 

Tabulation of data 

4 

16 

20 

30 

10 

Analysis 

60 

60 

60 

60 

30 

Hours/no, of sq 

114.5/ 

158.5/ 

260/ 

311/ 

131/ 

mile samples 

36 

144 

150 

250 

344 

Hours/squarc 

3.2 

1.1 

1.7 

1.2 

0.4 


mile 


RECOMMENDATIONS AND APPLICATIONS 

This investigation has addressed two primary areas of application for remote sensing techniques: 
1) the mapping and quantification of land use, and 2) the relationship of land use to stream sedimen- 
tation. 

The application for land use mapping based on aerial or satellite imagery has been well docu- 
mented in the literature. Conflicting demands for land utilization resulting from normal economic 
growth together with many legislative requirements indicate clearly that more data on land use 
will have to be generated and continually updated. The evaluation of a wide variety of imagery 
types in this study ted to the conclusion that high altitude color infrared imagery was the most 
universally acceptable for the mapping and measurement of land use types. Based on this Imagery 
a rapid and cost effective approach to the mapping and measurement of land use categories was 
demonstrated. The basic mapping approach developed in this investigation Ts now a part of the 
Corps of Engineers Handbook on Remote Sensing of the Environment (Haugen and Splett 1 974). 
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The comparison of taiui use lu stream sedimentation demonstrated that a relationship could be 
shown based primarily upon data derived by remote sensing means. This typo of data is necessary 
for the prediction of impact on streams and harbors of land use changes within a watershed. Thus, 
remote sensing techniques can solve a major problem in defining land usc>siream sedimentation 
relationships the sparsencss of published water quality data, especially for streams away from 
major population centers. 

Tire collection of data on water quality, especially sediment load, could be accomplished quickly 
and with minimum expense if a data bank were developed which would include incidcnt/rcficcted 
radiation measurements with each water sample made. With an adequate sampling over a regional 
area and throughout a representative range of national conditions, monitoring of stream sediment 
loads by remote sensing methods could be instituted as a practical and efneient technique. 
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Tabic AL Precipitation and temperature data far East Twin River Watershed* 
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APPENDIX A. PRECIPITATION AND TEMPERATURE DATA 
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Table AIL PrccipUalion and temperature data for Manitowoc River Watershed. 
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